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ABSTRACT 

We have carried out an extensive survey of magnetic field strengths toward 
dark cloud cores in order to test models of star formation: ambipolar-diffusion 
driven or turbulence driven. The survey involved ~ 500 hours of observing 
with the Arecibo telescope in order to make sensitive OH Zeeman observations 
toward 34 dark cloud cores. Nine new probable detections were achieved at the 
2.5-sigma level; the certainty of the detections varies from solid to marginal, so we 
discuss each probable detection separately. However, our analysis includes all the 
measurements and does not depend on whether each position has a detection or 
just a sensitive measurement. Rather, the analysis establishes mean (or median) 
values over the set of observed cores for relevant astrophysical quantities. The 
results are that the mass-to-flux ratio is supercritical by ~ 2, and that the ratio 
of turbulent to magnetic energies is also ~ 2. These results are compatible with 
both models of star formation. However, these OH Zeeman observations do 
establish for the first time on a statistically sound basis the energetic importance 
of magnetic fields in dark cloud cores at densities of order 10 3-4 cm -3 , and they 
lay the foundation for further observations that could provide a more definitive 
test. 

Subject headings: Stars: formation — ISM: magnetic fields — polarization 
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Introduction 



It has become increasingly clear that cosmic magnetic fields are pervasive, ubiquitous, 
and likely important in t he properties and evolutio n of almost everything in the Universe, 
from planets to quasars (jwielebinski fc Beck 2005 ). However, the role played by magnetic 
fields in star formation has remained uncertain. One extreme of the theories of star formation 
is that magnetic fields control the formation and evolution of the molecular clouds from 
which stars form, including the formation of cores and their gravitational collap se to form 



protos t ars. Detailed t heore tical work has been carried o ut by a number of groups. IShu et al. 



(119991 ) , iMcKed (119991 ) , and iMouschovas fc Cioleki (Il999l ) have reviewed and summarized the 
state of this theory. The fundamental principle is that clouds are formed with subcritical 
masses, M < M^,, where M^, = &/2tt\ /G, $ is the magnetic flux, and the expression for 
M$ is from lNakano fc Nakamura Jl978h . The magnetic field is frozen only into the ionized 
gas and dust; neutral gas and dust contract gravitationally through the field and the ions, 
accumulating mass (but not flux) in the cloud cores. This process is known as ambipolar 
diffusion. When the core mass reaches and exceeds M$, the core becomes supercritical 
(M > M$), collapses, and forms stars. The magnetic flux mostly remains behind in the 
envelope. Because the ambipolar diffusion time scale for the formation of cores is fairly long 
(> 10 7 yr), molecular clouds would have long lifetimes, althoug h observable "starless" core 
lifetimes would be short, ~ 10 6 yr (iTassis fc Mouschoviasll2004l ). 



This ambipolar-diffusion driven model for star formation was the standard for many 
years, but doubt of its validity w as raised by efforts to determine molecular cloud and core 
lifetimes. iBeichman et al.l (119861 ) used the number of molecular cores detected at various 
densities to determine typical statistical timescales for evolutionary sequences. Then the 
ratio of the number of molecular cores with T Tauri stars to the number of starless cores, 
together with typical T Tauri star lifetimes, led to an estimate of the li fetimes of star l ess co res. 
The result was a few times 10 6 yr. T his early estimate was refined bylLee fe Mversl (119991 ) to 
~ 0.3 — 1.6 x 10 6 yr. Later workers (IKirk et al.ll2005l ; iKandori et al.ll2005l ) inferred similarly 
short core lifetimes. This led to a new theory - that molecular clouds are intermittent 
phenomena, with short (~ 10 6 yr) lifetimes. In this theory clouds form at the intersection 
of turbulent supersonic flows in the interstellar medium. Generally, clouds do not become 
gravitationally bound, and they dissipate. Star formation occurs only in the smal l fraction of 
the molecular gas that is sufficiently dense to be self-gravitating (lElmegreenll2000l ) . Magnetic 
fields may be present in this theory, but they are too weak to be energetically important. The 
role of turbulence i n the e nergetics of the interstellar medium has been a very active area. 
Elmegreen fc Scalo ( 2004 ) have written an excellent review of interstellar turbulence, and 
MacLow fc Klessenl ( 120041 ) have extensively reviewed arguments that supersonic turbulence 
controls star formation. 
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Although the cloud lifetime estimates may seem to have doomed the strong magnetic 
field, quasi-static picture, that is not the case. First, cloud lifetimes seem to be longer 
than the free-fall time by a factor 2-5. Although supersonic turbulence may in principle 
provide support against collapse and lengthen cloud lifetimes, simulati ons have shown that 
such supersonic turb ulence will damp on roughly a free-fall time scale ( jMacLow et al.lll998 



Ostriker et al.lll999l). Something else, perhaps magnetic pressure, is slowing the collapse. 



Moreover, iTassis fc Mouschoviad (12004 ) have argued that the long cloud lifetimes that are 
quoted for the strong magnetic field model are for the entire lifetime from molecular cloud 
formation at density ~ 10 2 cm -3 to protostar formation, and that most of that time is spent 
in increasing densities to the values where cores can be identified observationally, ~ 10 4 
cm -3 . Once dense cores form, they argue that cores are magnetically critical to slightly 
supercritical, and that the core lifetimes in this model are 
observations. 



10 6 yr, in agreement with 



Hence, the evidence for cloud and core lifetimes does not provide a conclusive test. 
Because core lifetimes appear to be several times longer than free-fall times, there is evidence 
for a support mechanism in cores. However, this does not prove that magnetic fields provide 
that support. Moreover, the lifetime numbers are somewhat uncertain, due to the statistical 
nature of the arguments and the chain of reasoning from T Tauri star lifetimes to ages of 
cores at different densities. Finally, the core lifetime estimates do not address directly the 
timescale of the formation of dense cores from much lower density gas. 

The direct approach to resolving this uncertainly in the process by which stars form 
is to measure magnetic field strengths in molecular clouds in order to see whether they are 
weak (supercritical) or strong (subcritical). The crucial parameter is the ratio of the mass 
to the magnetic flux, M/$, which is of course closely related to M$. If M/$ is observed to 
be supercritical, particularly at lower densities, the magnetic support model is not viable. 
On the other hand, if it is observed to be subcritical, magnetic fields would be too strong 
for the intermittent, turbulent theory to hold. The M/$ parameter provides in principle a 
direct, simple, and definitive test to discriminate between the two theories of star formation. 



Of course, much work has been done to measure M/$ in molecular clouds. ICrutcher 



(119991 ) summarized the data available to that time, and con cluded that M/$ w as approx- 
imately critical to slightly supercritical in molecular clouds. iBourke et al.l (120011 ) extended 
the OH Zeeman work to the southern hemisphere with the Parkes telescope and added one 
and perh a ps two new detecti ons; they f ound results for Mj $ that essentially agreed with the 
Crutcherl (119991 ) conclusion. ICrutcherl (120071 ) has updated the discussion. However, the ex- 
tant observations have a major deficiency - a very small number of measurements of magnetic 
field strengths have been made in dark cloud cores - the sites of low-mass star formation. 
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Goodman et a 



Crutcher et al 



.1 (|1989f) detected the OH Zeeman effect toward Bl with the Arecibo telescope. 
( 19931 ) carried out a major and very sensitive OH Zeeman survey with the 
NRAO 43-m telescope toward 12 dark cloud cores, but they were only able to confirm the 
Bl detection and to obtain a possible detection toward p Oph. The large 18' beam of the 
43-m telescope meant that cores were not isolated; rather, the beam was filled by the cloud 
envelopes surrounding the cores. The 3' beam of the Arecibo telescope used for the survey 
reported here is well matched to the sizes of dark cloud cores a t distances of a few 100 pc . 
And in the early stages of the OH Zeeman survey reported here, I Crutcher fc Trolandl (120001 ) 
detected the OH Zeeman effect toward L1544. Hence, there are almost no measurements of 
magnetic field strengths in dark cloud cores. 

In this paper we report the full results of an extensive (~ 500 hours of telescope time) 
survey of the OH Zeeman effect toward dark cloud cores with the Arecibo telescope. In 
§2 we describe how the target cores were selected and the details of the observations. In 
§3 we present the results, including the inferred column densities, line-of-sight magnetic 
field strengths, and mass-to-flux ratios. In §4 we discuss these results, and in §5 present 
conclusions. 



2. Observations 



2.1. The Zeeman Effect 



Crutcher! (120071 ) reviewed the various techniques and results for studying magnetic fields 
in molecular clouds. Of these techniques, the Zeeman effect provides the only direct method 
for measuring magnetic field strengths in molecular clouds. In general, only those species 
with an unpaired electron will have a strong Zeeman splitting. This has limited detections 
to the the 21-cm line of H I, the 18- cm, 6-cm, 5-cm, and 2-cm A-doublet lines of OH, and 
the 3-mm N=1^0 lines of CN. The sole exception is the 1.3-cm H 2 maser line, due to very 
strong line strengths and strong fields in H2O maser regions. 

Except for some OH masers, the Zeeman sp litting is a small frac tion of the line width, 
and only the Stokes V spectra can be detected (ICrutcher et al.lll993l ); these spectra reveal 
the sign (i.e., direction) and magnitude of the line-of-sight component B/ os . By least-squares 
fitting the frequency derivative of the Stokes parameter I(v) spe ctrum dl{y)/dv to the 
observed V(u) spectrum, Bi os may be inferred (ICrutcher et al.lll993l ). 
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2.2. Selection of Targets 



Targets were selected from among molecular cores known to exist within the Arecibo 
declination range. Most targets fall into the Galactic anti-center direction (RA « 03 h - 07 h ); 
a few lie in the Galactic cen ter direction (RA 19 h - 21 h ). To identify suitable cores, we first 
examined CO survey data (IDame et al.ll200ll ) to identify nearby molecular cloud complexes 
observable from Arecibo. For each of these complexes, we then consulted detailed CO and 
NH3 maps of molecular cores to identify specific t argets. For example, in the anti-c enter 
reg ion, the L1457 cloud has been m apped in CO by IZimmermann &: Ungerechtd (119901 ) and 
by Moriartv-Schieven et al.l (119971). The Perseu s cloud has been surveyed for NH 3 cores 
by lLadd et al.l (119941 ) . lOnishi et al.l (119961 . 1 199a ) have identified and described a complete 
sample of 40 cores in th e TMC base d upon C 18 observations. The Mon OBI region has 
been s urveyed in CO bylOliyer et al.l (119961 ). And the Rosette Nebula cloud was studied in 
CO bv lSchneider et all Jl998h . 



Prior to Zeeman observations, we conducted a short survey of OH line strengths to- 
ward the principal cores in each cloud complex. We found that most cores in the Galactic 
center direction have relatively weak OH emission lines, owing, very likely, to the Galactic 
background continuum emission that is comparable to the excitation temperatures of the 
OH transitions. The final list of targets consisted of molecular cores for which the OH lines 
were relatively strong (since sensitivity to Bi os oc line strength); promising cores had to be 
deleted from the list due to lack of sufficient observing time. We did not want to spend very 
long times observing weak lines, for t hat woul d mean our survey would have few clouds. A 
sensitivity calculator such as that of iTrolandl (119901 ) allows one to estimate the observing 
time required to reach a given sensitivity. The actual sensitivity achieved varied due to vari- 
ations in line strength, line width, and available telescope time. Most targets are in nearby 
molecular clouds (e.g Perseus and Taurus) for which the Arecibo beam samples a relatively 
small linear scale (e.g. ~ 0.2 pc for a cloud at 200 pc). A few targets (e.g. the Rosette 
Nebula cloud) are more distant, but they add diversity to the sample since they are part of 
massive star-forming regions. 



2.3. Observations 



Simultaneous observations of the 1665 and 1667 MHz OH lines were carried out in 



the manner described by ICrutcher fc Trolandl (l2000l ). Very briefly, we used the single-pixel 
L-band wide receiver with native linear polarizations. A hybrid circuit immediately after 
the two HEMT receivers added ±90° of phase shift to the two linearly-polarized outputs to 
convert them into orthogonal circular polarizations. Typical system temperatures were 30- 
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40 K. Spectra were obtained with the "interim" correlator sampling 2048 spectral channels 
over a bandwidth of approximately 390 kHz, leading to a channel spacing of about 0.034 km 
s _1 . Spectra were Hanning smoothed to produce a final spectral resolution of two channels. 
Data analysis followed our usual procedures for Arecibo OH Zeeman observations. We 
constructed Stokes I and V profiles from the sum and difference, respectively, of the line 
profiles in orthogonal circular polarizations. Then we least-squares fitted to each V profile 
a function defined by the sum of a constant times the I profile plus another constant times 
the derivative of the I profile. The first constant is a measure of gain differences between 
the two circular polarizations; this constant was found to be negligible for all sources. The 
second constant is proportional to the line-of-sight magnetic field component Bi os . Finally, 
we computed a weighted average of Bi os as derived independently from the two OH emission 
lines. The full survey results are given in tables [T] and [2J Figure [T] shows the Stokes I and V 
spectra for the 9 positions with a Zeeman signal > 2.5cr, and figure [2] shows all of the inferred 
results for B; os from the Arecibo survey plotted against N(H2). Six of the nine positions in 
figure [H have solid Zeeman detections. In order of increasing marginality of detection for the 
remaining probable detections, there are B217-2, B5, and L1457Sn. See detailed discussion 
below and in the appendix. 

We carried out several tests to establish the reliability of the Arecibo system for Zeeman 
observations. For one, we established the sense of circular pola rization with two helical tes t 



antennas of known circular polarization sense, as described by lCrutcher fc Trolandl (120001 ) . 
We also observed the highly circularly polarized masers in W49 and compared the Stokes 
V profile with those in the previous literature. Over the duration of the project, we often 
observed the weak but highly circularly polarized OH maser in S247 to establish that no 
significant change in circula r polarization response had occurred. Finally, as also described by 



Crutcher fc Trolandl (120001 ). we estimated beam squint via continuum polarization mapping 
observations of an unresolved continuum source. These data establish that the difference in 
pointing between the right and left circularly polarized beams (i.e. the beam squint) was of 
order 1", less than 1% of the approximate 3' FWHM beamwidth of the Arecibo telescope at 
1666 MHz. 

In addition to the Zeeman observations of molecular cores, we mapped OH line strengths 
in the vicinities of the Zeeman positions. These mapping observations serve two purposes. 
For one, they allow us to determine the sizes of the cores on the plane of the sky, from which 
core volume densities and masses may be estimated (tables 1 & 2). Our crude OH maps 
were supplemented by C 18 maps from the literature. We used volume density n(H 2 ) = 
N(H 2 )/2r and mass M = 7rr 2 A^(if 2 )2.8m^ to estimate these parameters, where r is the 
cloud radius and the factor 2.8 includes a 10% He abundance and m# is the mass of an H 
atom. Also, the mapping observations permit us to establish upper limits for each core upon 
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the instrumental effects of beam squint. For this latter purpose, we used mapping positions 
offset approximately 3' (one beamwidth) north, south, east and west from each Zeeman 
position. From 1667 MHz profiles at these offset positions, we estimated the magnitude and 
direction of the linear velocity gradient in the OH emission at the Zeeman position. Then we 
computed the expected instrumental magnetic field from a beam squint of 1" if the position 
angle of the beam squint exactly matches the position angle of the velocity gradient on the 
sky. We feel this procedure offers a realistic upper limit to the instru mental Zeema n effect, 
even if the actual beam squint is, in practice, sli ghtly higher than 1". (IHeilesI (119991 ) reports 
beam squints of typically 1.3" at 1420.4 MHz: ICrutcher fc Trolandl (boooh report a 1666 
MHz beam squint of 1.2" ± 0.2"). The actual velocity gradient is unlikely to be aligned with 
the telescop e beam squint, esp ecially since the beam squint position angle is nearly fixed 
in azimuth (IHeiles et al.ll200ll ). Therefore, the beam squint will rotate on the sky as the 
parallactic angle of the source changes during a Zeeman source observation. As a result, the 
effective beam squint over a several hour Zeeman observation will be significantly less than 
the actual beam squint. 

Estimates of upper limits to instrumental effects (see above) suggest that the magnetic 
field detections reported here are reliable. The upper limits to instrumental effects have an 
average value of 2.4 /zG for all of our Zeeman positions. The highest instrumental field is 
5.4 fiG at position Ros4, for which cr(B) is 9.8 fiG. Of the sources listed in table 2, nine 
have magnetic fields with S/N > 2.5. These sources are likely to be Zeeman detections 
on statistical grounds alone (see §3). In none of these sources is the upper limit to the 
instrumental magnetic field more than 25% of the derived field value (table 2, column 4), 
for most sources it is less. We conclude that the magnetic field detections (i.e., S/N > 2.5) 
are quite unlikely to result from instrumental effects associated with beam squint. More- 
over, non-detections but sensiti ye measurements a re unlikely to be significantly affected by 
instrumental effects. Note that IHeiles et all (120011 ) report that beam squint is the principal 
contributor to instrumental circular polarization with the Arecibo L-band wide receiver. 



3. Results 

3.1. Measurements versus Detections of Bi os 

Although the astrophysical analysis we will report on below uses each of the 34 measure- 
ments of Bi os without directly considering whether each measurement is a detection or not, 
it is of interest to consider which of the 34 measurements should be regarded as detections. 
Since detections are not produced by instrumental effects (§2), there are two statistical cri- 
teria and a third subjective criterion: (1) that \Bi os \/ ctb 1os be greater than an appropriate 
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number, (2) that the 1665 and 1667 MHz lines yield the same result for Bi os within the 
measurement uncertainties, and (3) that a plot of Stokes V looks consistent with detection 
of the Zeeman effect. 

For criterion (1), we set \Bi os \/ <jb 1os > 2.5, the level at which there are probably no 
false detections for a sample size of 34. The normal probability function says that for a 
normal error distribution, the fraction of the measurements that should be 2.5cr or more 
from the "real" value is 0.0124. In our case, we have 34 positions at which we measured 
Bi os ; 0.0124 x 34 = 0.4. So a criterion of 2.5<j with a sample of 34 means we would claim 0.4 
false detections. Based on this criterion alone, there are 9 detections; the probability is that 
there are no false detections, but it is possible (less than 50% probability) that there is 1. 

For criterion (2), we set AB/a AB < 1.9, where AB = |£ tos (1665) - B los {1667)\. The 
normal probability function says that the fraction of the AB that should be 1.9cr or more 
from zero is 0.057, and with 9 possible detections, this criterion would yield 0.057 x 9 = 0.5 
possible false positives. The largest AB/aAB in table[2]is 1.7. Hence, all 9 possible detections 
meet this consistency criterion. 

Finally, there is the subjective criterion (3). Figure [1] shows Stokes I and V plots of 
each of the 9 probable detections. For these plots we have combined the 1665 and 1667 MHz 
results, weighted by the inverse square uncertainties in Bi os for each transition. The Stokes I 
spectra are the weighted sum of Ta(1665) and 7^(1667); these are then the Stokes I spectra 
that would have been observed if there were a single OH line rather than two. The Stokes 
V spectra are the observed Stokes V spectra for such a single line. This combination gives 
the appropriate spectra for judging the significance of each probable detection. 

Although 9 positions pass our criteria for detection, clearly some are more solid than 
others. We discuss each of the 9 possible detections in the appendix with respect to the 
likelihood that each is or is not a detection. Again, however, note that our analysis does not 
depend on whether a particular position does or does not have a detected Zeeman signal. 



3.2. Estimating B tota i and M/$ from B, 



log 



Following [Crutcheij (119991 ). we define 

A = (M/$) ofeOTd /(M/$) 

critical > 



(1) 



where (M/Q) b serve d is the observed mass to flux ratio inferred from the r atio of iVfi^) to B 
and ( Mj <3>) critical is the theoretically determined critical mass to flux ratio (INakano &: Nakamura 



19781 ). 



- 9- 



Then 



A = 7.6 x 1Q- 21 N(H 2 )/B. 



OS ■ 



(2) 



where N(H 2 ) is the column density of H 2 in cm 2 and Bi os is the line-of-sight magnetic field 
strength in /zG. 

It must be kept in mind that all of the B/ os results are lower limits to the total magnetic 
field strength. It is possible to correct statistically for the fact that only one component 
of B is measured, i.e., B^ os = |B| cos 6*. For a large number of clouds with the same total 
field strength for which the angle between B and the observed line of sight is randomly 
distributed, 



~T5 _ -/o 

Dins — 



tt/2 



B I cos sin 9d9 



f* /2 sin 9d9 



I|B|. 

2 1 1 



(3) 



If B is strong, clouds w ill have a disk morphology with B along the minor axis (cf, 
Mouschovas fc Cioleki (119991 )). To properly measure A, one needs B an d N along a flux 
tube, i.e., parallel to the minor axis. Then, as noted by ICrutcherl ( 119991 ). the path length 
through a disk will be too long by 1/ cos# and iV will be overestimated, while |B| will be 
underestimated by cos 9. Statistically, 



— — C /2 (M/<S>) obs cos 2 9sm9d9 1 , s 
J sm Odd * 



(4) 



3.3. Other Physical Parameters 



the OH 


data. 


( Crutcher 


1979) 



a x T x AV x 10 



ii 



cm -2 K" 1 km s" 1 , where T is the peak line antenna temperature, AV is the full-width-at- 
half-maximum intensity of the line, a = 8.49 for the 1665 MHz line and 4.71 for the 1667 MHz 
line. The coefficients include the beam efficiency of the telescope, t]b ~ 0.5. The N(OH) in 



table 1 are the average results from the two lines. Then N(H 2 ) = N(OH)/8x 10 8 (j Crutcher 



19791 ). 



A comparison of the total masses estimated from N(OH) and r (which we designate 
as Moh) with the virial masses (M^riai) may be instructive. The means for all 34 cores 
are M Q h = 58 M and M viria[ = 74 M ; the respective medians are M Q h = 16 M and 
M V i r i a i = 29 Mq. Considering the mean values, the ratio M OH /M virial = rN(OH)/AV 2 X = 
0.78, where X is the OH/H abundance ratio. The mean mass ratio would be 1 if any one 
of the following were true: r larger by 1.3, N(OH) larger by 1.3, AV smaller by 0.88, or 
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OH/H = 6 x 10~ 8 . Any or some combination of all of these is possible. The r we use may 
not be the relevant r, both because of the way OH samples the gas and the geometrical 
assumption (§2.3) made in going from iV and r to Moh- N(OH) could be larger if the lines 
were slightly saturated rather than being optically thin as assumed. The relevant AV for the 
virial mass may be different from that given by OH, again because of how OH s a mples the H 2 



gas distribution. And the value of OH/H may be different from the ICrutcherl (119791 ) result. 
Moreover, the virial mass calculation assumes virial equilibrium between gravity and kinetic 
motions, which may not be correct. The cores may not be in equilibrium, and the support 
provided by the magnetic field is not included. However, rather than attempting to discuss 
these points, we prefer to emphasize that the agreement in the two methods of determining 
mass is quite good; our estimates for astrophysical quantities such as the mass-to-flux ratio 
and the Alfvenic Mach number are therefore unlikely to be dominated by systematic errors 
in column and volume densities. 

Because we can only measure one component of the magnetic vector B, it will be 
necessary to consider the mean or median values of the mass-to-flux ratio, ratio of turbulent- 
to-magnetic energy, and the Alfvenic Mach number for the ensemble of dark cloud cores we 
observed. The values for the parameters that go into calculating these quantities are the 
following mean and median values, respectively: Av = 0.71, 0.71 km s -1 , Bi os = 8.2, 5.8 pG, 
r = 0.31,0.29 pc, N 21 {H 2 ) = 4.5,4.0 cm" 2 , and n(H 2 ) = 3200, 1800 cm" 3 . These values will 
be used in the discussion below. 



4. Discussion 



An important motivation for this Zeeman effect survey was to better define the key 
processes responsible for star formation, in particular, the role of the magnetic field. The role 
of the magnetic field in star formation, in turn, depends upon the ratio of magnetic energy 
in star forming (i.e. self-gravitating) clouds to other relevant energies. These energies are 
the gravitational energy of the cloud and the energy of internal motions. Since star forming 
clouds are, in general, not rotating nor collapsing, the energy of internal motions is presumed 
to reside in macroscopic motions usually referred to as turbulence and measured by the line 
widths of molecular spectral lines. The parameter A, described in §3.2, is a measure of the 
ratio of gravitational to magnetic energies in the cloud. The ratio of turbulent to magnetic 
energies is given by (3 tU rb = 47rF pa 2 / B 2 , where p is the density, a is the 1-D turbulent 
velocity dispersion, B is the total magnetic field strength, and F=3 for 3D turbulence. Then 
(3 turb = 0.32n(H 2 )[AV 2 - 0.027]/3Bf M , where n(H 2 ) is in cm" 3 , AV is the FWHM in km 
s" 1 , the line-of-sight field strength Bi os is in pG, and 0.027 removes the thermal contribution 
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to the observed line widths for an assumed kinetic temperature of 10 K. 

Values for these key ratios (A and I3 tur b) are not meaningful for individual clouds since 
we cannot measure the total field strength, only Bi os . It is for just this reason that a survey 
of many clouds is necessary to adequately characterize the role of the magnetic field in these 
objects. Therefore, we use the data from the present survey to estimate mean and median 
values for the magnetic field strength and for A and f3 tU rb over the ensemble of clouds observed 
for this survey. 

In order to proceed with the estimation of average quantities, we make the assumption 
that the angle between the line of sight and the magnetic field is randomly distributed 
among our 34 positions. For this assumpt ion, the probability den sity function for Bi os is 



flat between and the total strength |B| (IHeiles fc Trolandl 120051 ) . That is, the observed 



Bi os are assumed to differ from a single total field strength for the sample of cores only 
by geometrical projection. We can then calculate the weighted mean Bi os and therefore A, 
where the weighting is by the inverse square uncertainties in Bi os ; the uncertainty in A is 
dominated by the measurement uncertainties in Bi os rather than in N(OH). The result 
is Bios = 8.2 ± 2.2 /iG and A ~ 4.2; this is shown as a solid line in figure [2j If median 
values are used instead, we find the median value of the mass-to-flux ratio A1/2 ~ 5.2. 
All measurements, detections and non-detections, are included in the weighted mean and 
median. Thus, it makes no difference for these calculations whether a given position has a 
detected Bi os ; what matters is each measured value and its uncertainty. Our result for the 
total mean field strength (equation [3]) is then |B| = 16.4 /zG. The above value for A is then 
systematically too large due to geometrical effects. The minimum correction is a factor of 
1/2 due to the correction for only measuring B[ os . For a disk morphology, the correction is 
1/3 (equation |4]) . Hence, the geometry corrected A c ~ 1.4 — 2.1; Ai/2, c = 1-7 — 2.6. Finally, 
f3 turb = 2.4 and (3 tU rb,i/2 = 2.7. The related Alfvenic Mach numbers are Ma = 1-4 and 
M A1/2 = 1.5. 

Also plotted in figure [2] as a dotted line is the critical mass-to-flux ratio. Although 
we measure only the line-of-sight component of the magnetic field, one might expect the 
magnetic field to point approximately along the line of sight in a few cases. If there were a 
subcritical core in our sample and the field pointed along the line of sight, its plotted value 
would lie above this critical line. 

What conclusion can we draw from the estimated value of A c or Ai/2, c ? The cores 
of molecular clouds (sampled at a typical density of a few times 10 3 cm -3 ) appear to be 
slightly supercritical by about a factor of 2. Therefore, on average the gravitational energies 
somewhat exceed magnetic energies in these clouds, although the uncertainties in the results 
do not rule out a critical mass-to-flux ratio. Had the estimated value of A c been significantly 
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higher than 2, we could have ruled out the ambipolar diffusion model of star formation 
since magnetic fields, in this circumstance, would be energetically insignificant. Likewise, 
had the estimated value of A c been less than unity, we could have ruled out the turbulent- 
driven model owing to the strong influence of magnetic fields on molecular cores. As it 
is, the estimated value of A c is consistent with both extreme-case models. The ambipolar 
diffusion model predicts that cores are formed in subcritical clouds, but by the time cores 
with n(H 2 ) ~ 10 4 cm -3 have formed by the action of ambipolar diffusion, they are critical 
to slightly supercritical. The turbulence driven (i.e. weak magnetic field) model forms cores 
by turbulent compression. Although these models have the envelopes and regions between 
cores highly supercritical, many cores may be only slightly supercritical. Hence, although 
the observations have shown that magnetic fields are sufficiently strong that they cannot 
be ignored, the very hard-won observational results cannot rule out either model of star 
formation. 

What conclusion can we draw from j3 turb = 2.4? Obviously, this value suggests that 
turbulent energy exceeds magnetic energy, just as gravitational energy exceeds magnetic 
energy (although not overwhelmingly in either case). So the cloud cores sampled in this 
project, on average, are slightly out of equipartition between turbulent and magnetic energies. 
Note that l3 tUT b is related to the Alfvenic Mach number by the relation M\ « f3 tU rb- Therefore, 
we infer a mean value M a ~ 1-6. That is, internal motions are on average mildly super- 
Alfvenic in these cores. If turbulent energy exceeds magnetic energy by a factor of a few, (and 
the cores are not collapsing), then the cores are in approximate virial equilibrium between 
internal motions and gravitation, with magnetic support of lesser significance. The average 
ratio of mass to virial mass in table 2 is 0.8, a value close to unity and consistent with this 
conclusion. 

A concern that sometimes arises regarding Zeeman effect observations (and the conclu- 
sions drawn there from) involves tangling of the magnetic field on scales smaller than the 
beam. Field tangling, of course, can reduce the beam averaged value of Bi os ; a field reversal 
within the beam can conceivably reduce the value to 0. Although the Arecibo beam is small 
compared to beam sizes used for some other Zeeman effect projects, field tangling can still 
exist, in principle, on scales smaller than any beam. Nonetheless, we do not believe that 
field tangling on scales smaller than the Arecibo beam has significantly affected this project, 
nor the conclusions we draw from it. For one, higher spatial resolution studies of magnetic 
fields in molecular cores show little sign of field tangling. These studies rely upon linear 
polarizat ion of dust emission to map th e morphology of the field in the plane of the sky. For 



example, IWard-Thompson et al.l (120001 ) mapped linear polarization in L1544 (also included 



in this project), L183 and L43 at 14" spatial re s olutio n. They find a high degree of order 



in the field. Likewise, iGirart. Rao, fc Marrond (120071 ) find a high degree of order in the 
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field of NGC 1333 IRAS 4A at 1.5" spatial resolution. Another consideration regarding field 
tangling is the fact that the value of A is related to the net magnetic flux through a cloud. 
To first approximation, at least, the beam-averaged value of Bi os is a measure of just that 
quantity (for a single field component), regardless of any small-scale field irregularities that 
may exist on scales smaller than the beam. Finally, the data from this project suggest that 
the observed field strengths are sufficient to be energetically important (although not domi- 
nant) in the molecular cores. Therefore, the observed field strengths are sufficient to impose 
a reasonable degree of order on the field. If the observed field strengths are significantly 
lower than the actual field strengths, owing to tangling, then the actual magnetic field must 
be even more energetically important. In such a case, field tangling is even less probable. 
That is, the hypothesis that field strengths in molecular cores are higher than implied by 
this project, owing to tangling, leads to a conclusion that contradicts the hypothesis itself. 



5. Conclusions and future directions 

This project had the potential to distinguish between the two extreme models, ambipo- 
lar diffusion versus turbulence driven star formation. However, the result for A c lies in a 
range compatible with both models since the molecular cores surveyed appear to be slightly 
magnetically supercritical. That is, gravitational energy appears to dominate magnetic en- 
ergy by a small factor. Also, the estimated value of f3 turb ~ 2.4 is such that turbulent energy 
in the cores appears to dominate magnetic energy, again, by a rather modest factor. To 
zeroth order, at least, gravitational, magnetic and turbulent energies in these cores appear 
to be comparable. Once again, nature appears to have held her cards regarding magnetic 
effects on star formation close to her chest! At the same time, the OH Zeeman observations 
reported here are of considerable value since they do establish for the first time on a statisti- 
cally sound basis the energetic importance of magnetic fields in dark cloud cores at densities 
of order 10 3-4 cm -3 . 

What future observations might help to distinguish between the two extreme models 
of star formation? Further Arecibo observations like the ones reported here are unlikely 
to settle the issue, especially given the observing time that would be required to improve 
significantly these results. Instead, there are at least two types of Zeeman effect observations 
that have the potential to do so. One type of observation targets magnetic field strengths 
in the envelopes of clouds, not in their cores. The ambipolar diffusion-driven model of star 
formation predicts that these envelope or inter-core regions must be magnetically subcritical 
(A < 1). The turbulence-driven model predicts that the same regions must be magnetically 
supercritical. Zeeman effect measurements of OH emission lines in the inter-core regions 
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of molecular cores would be very difficult owing to the weakness of the OH emission there. 
However, OH absorption lines against extra-galactic continuum sources have the potential to 
probe the field along random lines of sight through the clouds. Only the Arecibo telescope 
has the ability to perform this experiment since its large collecting area ensures an adequate 
number of background continuum sources to make the observations meaningful. Another 
approach to distinguishing between the two models of star formation is to measure differential 
mass-to-flux ratios, that is, to measure A in cloud cores relative to A in the envelopes of the 
same clouds. The ambipolar diffusion model predicts that A is greater in the cores than 
in the envelopes, owing to the ambipolar diffusion process itself (see §1). The turbul ence 
drive n model predicts the opposite, namely, that A decreases with increasing density (IDib 



20061 ). Therefore, it should be possible to distinguish between the two models by measuring 



A in the core and envelope of the same cloud. The OH Zeeman effect offers an opportunity 
to make such measurements. This approach has the advantage of eliminating geometrical 
effects associated with the measurement of Bi os , since, to first approximation, the angle 
between the field and the line-of sight should be the same in both core and envelope of the 
same cloud. We are currently pursuing both approaches (OH absorption and differential 
A measurements) in an attempt to better distinguish between the alternate models of star 
formation. 



We thank the Arecibo Observatory for the generous allocation of observing time, which 
allowed this project to be carried out, and the observatory staff for help in making the 
observations successful. This research was partially supported by NSF grants AST 0205810, 
0307642, and 0606822. 



A. Discussion of possible detected Zeeman signals 

L1457S: The line has a single strong peak and an extended wing to positive velocities. 
The Stokes V spectrum shows that the Zeeman signal is clearly detected from the strong 
peak. 

L1457Sn: This short integration was obtained 3' north of the L1457N position to look 
for spatial variation in Bi os . The line profile is qualitatively similar to that toward L1457S. 
Although the fitting process reported a marginally stronger field at this position than at the 
L1457S position, the noise level in the Stokes V data was too high to see the Zeeman signal. 
In order to enhance its visibility, the observed Stokes V spectrum was boxcar smoothed 
by 5 channels; that is the spectrum that is shown. The displayed fit line has also boxcar 
smoothed. Since the integration time at this position is only about l/9th that toward L1457S, 
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this smoothing still has a higher channel-to-channel noise level than the unsmoothed L1457S 
Stokes V spectrum. Yet there is still not a clear Zeeman signal in the observed Stokes V 
spectrum. The fit seems to be responding to the 8 of 9 channels that are negative to the 
negative velocity side of the line peak and the 2 positive channels to the high velocity side. 
Nonetheless, we regard L1457Sn as having failed the subjective criterion (3) test. On the 
other hand, this position is only 1 FWHM beam away from the L1457S position; the Bi os at 
the two positions have the same sign (field direction) and consistent magnitude. In spite of 
failing criterion (3), we still regard the measurement toward L1457Sn as being more likely a 
detection than not, although a rather marginal one. 

L1448CO: The line profile is clearly asymmetric, suggesting that there are multiple 
velocity components. The fit responds most strongly to the steeper slope on the lower 
velocity side of the line and the negative feature in Stokes V. The observed Stokes V appears 
to have a corresponding positive feature just to the positive velocity side of the peak, but 
the dl/du fit does not fit this feature well, due to the smaller slope of the observed Stokes 
I spectrum at these velocities that appears to be caused by a weaker additional velocity 
component at these velocities. The detected B\ os is confined to a narrow velocity component 
that has only its more negative side (apparently) unaffected by line component blending. 

L1448COe: Observing this position, that is 3' east of the L1448CO position, was again 
an attempt to look for small-scale structure in B\ os . Here, unlike L1456Sn, there is a clear 
detection of the Zeeman signal. The Stokes I line profile shape is somewhat different from 
that at the L1448CO position, being at slightly more positive peak velocity and lacking the 
steep slope on the more negative velocity side. The close agreement in Bi os between the 
two positions suggests that in spite of the difference in the Stokes I profile, there is little 
difference in Bi os . 

B5: B5 has an apparently simple, single line profile. The Zeeman signal in the Stokes V 
spectrum is not unambiguous, but is consistent with the marginal detection given by the fit. 
Although this source passes our three tests, the detection should be regarded as marginal or 
probable rather than certain. 

B217-2: The Stokes I profile is fairly simple, although there is wide, weak plateau 
emission that does not, however, interfere with the Zeeman signal, which appears clear and 
unambiguous. The less than 2a detection of the Zeeman effect in the 1665 MHz line could 
lead to this being considered a marginal detection, although it fully meets our three detection 
criteria. 

TMC1: There is clearly a complex profile, with at least three strong narrow components. 
Only the lowest velocity component has a clear Zeeman signal, at least on the low velocity 
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side. The higher velocity side is confused with the middle velocity component, and there 
is no Zeeman signal from the higher velocity line components. This does not necessarily 
mean there is only a very week field in these components; blending may obscure the Stokes 
V Zeeman signatures from these components. 

L1544: This starless core has a fairly simple line profile and an unambiguous Zeeman 
signal in the Stokes V spectrum. 

L Oril: This position clearly has at least two velocity components, with the unambigu- 
ous Zeeman signal coming from the stronger, more positive velocity one. 
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Table I. Dark Cloud Core OH Results 
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c Distance, pc 

d mean cloud radius to OH half-power point, pc 
e hours of integration time 
! Vlsr, km s _1 

s full- width at half-maximum, km s _1 
h H 2 column density, H 2 x 10 21 cm~ 2 
'H2 volume density, H2 cmT z 
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Table 2. Dark Cloud Core Zeeman Results 
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Fig. 1. — Arecibo Stokes I and V spectra for the 9 probable detection positions in the survey. 
Observed data are histogram plots; fits to Stokes V are the dark lines. These are weighted 
means of the 1665 and 1667 results. The Stokes V spectra have been scaled up and shifted 
by -0.3 K for display purposes. 




Fig. 2. — Results for B; os from the Arecibo dark cloud survey plotted against the H 2 column 
density (N 2 i = 10~ 21 N). The 9 probable detections (see text) are plotted as filled circles, 
while non-detections are plotted as open circles. Error bars are la. The solid line is the 
weighted mean value for the mass to flux ratio with respect to critical inferred from the 
Zeeman Bi os data with no geometrical correction; A 4.8±0.4. After geometrical corrections 
(see text), A c 2, or slightly supercritical. The dashed line is the critical mass to flux ratio. 



